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Abstract. For two-dimensional Euler equation on the torus, we prove that 
the L°°-norm of the vorticity gradient can grow as double exponential over 
arbitrary long but finite time provided that at time zero it is already sufficiently 
large. The method is based on the perturbative analysis around the singular 
stationary solution studied by Bahouri and Chemin in pQ. Our result on the 
growth of the vorticity gradient is equivalent to the statement that the operator 
of Euler evolution is linearly unbounded in Lipschitz norm for any time t > 0. 



1. Introduction and some upper bounds 

Consider the two-dimensional Euler equation for the vorticity 

6 = V6-tP, iP = V ± A- 1 9, 6(x,y,0) = 6 (x,y) (1) 

and 8 is 27r-periodic in both x and y (that is, the equation is considered on the torus 
T 2 ). We assume that 6*o has zero average over T 2 and then A^ 1 is well-defined since 
the Euler flow is area-preserving and the average of 9(-,t) is zero as well. Denote 
the operator of Euler evolution by £t, i.e. 

0(f) = £ t 9 

The global existence of the smooth solution for smooth initial data is well-known 
and is due to Wolibner [T3] (see also [TO])- The estimate on the possible growth 
of the Sobolev norms, however, is double exponential. We sketch the proof of this 
bound for H 2 - norm. The estimates for H s ,s > 2 can be obtained similarly. More 
general results on regularity can be found in [3]. Let 

hit) = 

Lemma 1.1. If 9 is the smooth solution of {Ip, then 

m < exp (a + 21o g + . 2 (0))exp(C7||^o||^)-l^ (2) 

Proof. Acting on ([T]) with Laplacian we get 

Af? = A9 x tp y + 2V9 X ■ V4' y - A9 y ip x - 2V9 y ■ Vip x 
Multiply by A9 and integrate over T 2 to get 

dt\\0\\ 2 H2 < ||if(VOHoo||0|| 2 H 2 (3) 
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where H(ip) denotes the Hessian of ip. The next inequality follows from the 
Littlewood-Payley decomposition (see [4], proposition 1.4 for more general result) 

UflXVOIloo < C(a)||0||oo(l + log+ \\6\\ H .) (4) 

for any a > 1. Notice that ||0||oo is invariant under the motion so combine ([3]) and 
Q to get ©. □ 

Remark 1. In the same way one can prove bounds for higher Sobolev norms, 
e-g-, 

logj 4 (i) < (1 + log + j 4 (0))exp(q|(?o|M) - 1 (5) 
The natural questions one can ask then are the following: first, how fast can the 
Sobolev norms grow in time and what is the mechanism that leads to their growth? 
Secondly, for fixed t, how does ||£(t)0 o ||i?» depend on ||^oIIh s when the last expres- 
sion grows to infinity? For example, given ||#ol|oo ~ 1; the right hand side in @ 
grows as a power function in ,72(0), the degree grows exponentially in t and is more 
than one for any t > 0. 

Instead of working with Sobolev norms, we will be studying the uniform norm 
of the vorticity gradient (or Lipschitz norm) as this norm is more natural for the 
method used in the proof. It allows the similar upper bound. We again give the 
sketch of the proof for completeness. 

Lemma 1.2. If 0q is smooth and ||#o||oo ~ 1 ; then 

||V£ t o ||oo < exp (C(l + log+ ||V0 o ||oo)e Ct ) (6) 

Proof. If ty(z,t) is area-preserving Euler diffeomorphism, then 

(£ t 9 )(z) = 6 (*- 1 (z,t)) 

On the other hand, V&(z,t) solves 

* = -ii(*,£), *(z,0) = z 

where u(z,t) = V ± A~" 1 0(-,t). For the Riesz transform we have a trivial estimate 

(A-^Hoo < 1 + U0IU1 + log+ HV^Iloo) (7) 

( Indeed, without loss of generality we can evaluate the integral at zero and assume 
0(0) = 0. Then, e.g., 



JBi(O) ISI JB s (0) J8<\z\<l ISl 

where <5 _1 = max{|| V0||oo, 2}. Apply now the Lagrange formula to the first term 
to get 0. 



so 

\u(w u t) - u(w 2 ,t)\ < \wx-w 2 \b, & = l + log+||V0(i)||oo 
Therefore, we have 

l/l<A f(t) = \*(z 2 ,t)-*(z u t)\ 2 , 

After integration 

|z 2 -2i|expf-C^ b(r)dr^ < \*(z 2) t)-*(zx,t)\ < \Z2-Zi\expfc b(r)dr 
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Since 



||V0(«,t)||oo = sup 



\Z2 ~ Z\\ 



we get inequality 




V Jo 



Taking log of the both parts and applying the Gronwall-Bellman, we get ([6]). 



In this paper, we will work only with large ||V#o||oo- For that case, we will show 
that, given arbitrarily large A, the estimate max te r 0;r i ||V0(-, i)||oo > A eT_1 || V^olloo 
can hold for some infinitely smooth initial data. This is far from showing that ([2]) 
or ^ are sharp however it already is equivalent to the statement that £( is linearly 
unbounded. The question of whether ||V#o||oo can be taken ~ 1 is left wide open, 
see discussion in the last section. 

Our results rigorously confirm the following observation: if the 2D incompressible 
inviscid fluid dynamics gets into a certain "instability mode" then the Sobolev 
norms can grow very fast in local time (i.e. counting from the time the "instability 
regime" was reached) . Can the Sobolev norms grow at all infinitely in time assuming 
that initially they are small? The answer to this question is yes, see [5] and [51IT21I51 
1141111] . The important questions of linear and nonlinear instabilities were addressed 
before (see, e.g., [3] and references there). In the recent paper [7J, it was proved 
that £ t is not uniformly continuous on the unit ball in Sobolev spaces. 

Remark 2. It must be mentioned here that 2D Euler allows rescaling which 
provides the tradeoff between the size of 6q and the speed of the process, i.e. if 
9{x, y, t) is a solution then fi9(x, y, fit) is also a solution for any \i > 0. However, in 
our construction we will always have ~ 1, Vp € [l,oo]. 

Remark 3. If one replaces A^ 1 in (TTJ) by A~ a with a > 1, then the growth of 
the vorticity gradient is at most exponential, e.g. 



Moreover, the lower exponential bound can hold for all times as long as 9$ is 
properly chosen (see [5]). 

2. The singular stationary solution and dynamics on the torus 

The following singular stationary solutions was studied before (see, e.g., [TJ|3] in 
the context of R 2 ). We consider the following function 



This is a steady state. Indeed, the function ipQ — A^ 1 ^ is odd with respect to 
each variable as can be verified on the Fourier side. That, in particular, implies 
that ipQ is zero on the coordinate axes so its gradient is orthogonal to them. This 
steady state, of course, is a weak solution, a vortex-patch steady state. Another 
consequence of tpo being odd is that the origin is a stationary point of the dynamics. 
By the Poisson summation formula, we have 



□ 



l|VSr J dolU < ||Vflo||ooexp(C||tf ||oo*) 



Sq(x, y) = sgn(x) ■ sgn(y), \x\ < tt, \y\ < tt 



C\n\z\+<j){z) 7 z~0 



neZ 2 ,ra/(0,0) 
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where <fi(z) is smooth and even. 

Therefore, around the origin we have 

VMx, y) ~ / / , {X ~S U ",~H u *9n(Si)sgn(Z 2 )dZidS2 + (0(y), 0(x)) 

Due to symmetry, it is sufficient to consider the domain D = {0 < x < y < 0,001}. 
Then, taking the integrals, we see that 

y) = (Mi,M2) = (V^Vo) (a;, y) = (8) 

= d(-J ln(y 2 +e 2 )^ + xn(x,2/),^ ln(x 2 + f)d£ + yr 2 (x, y)j 

= c 2 (-x\ogy + xO(l),y\ogy + yO(l)) if (x,y) e D 

The correction terms ru 2 ) are smooth. Without loss of generality we will later 
assume that c 2 = 1 in the last formula (so C\ — 0.5). That can always be achieved 
by time-rescaling. Notice also that the flow given by the vector-field fi is area- 
preserving. 

Thus, the dynamics of the point (a, f3) G D, a <C f3 is 

{CiPf < y(t) < (C 2 (3) e \ a{C l pr Et+1 < x(t) < a{C 2 p)- Et+ \ t e [0,t ], (9) 

where to is the time the trajectory leaves the domain D. These estimates therefore 
give a bound on t . The attraction to the origin, the stationary point, is double 
exponential along the vertical axis and the repulsion along the horizontal axis is 
also double exponential. 

3. The idea 

The idea of constructing the smooth initial data for a double exponential scenario 
is quite simple and roughly can be summarized as follows: given any T > 0, we will 
smooth out the singular steady state such that the dynamics is double exponential 
over [0,T] in a certain domain away from the coordinate axes. Then we will place 
a small but steep bump in the area of double exponential behavior and will let it 
evolve hoping that the vector field generated by this bump itself is not going to ruin 
the double exponential contraction in OY direction. The rest of the paper verifies 
that this indeed is the case. 



4. The Model Equation 

Consider the following system of ODE's 

x = fj, 1 (x,y) + u 1 (x,y,t), x(a,/3,0)=a 
y = fj, 2 (x, y) + v 2 {x,y,t), y(a,P,Q) = P 

Here we assume the following 

I f 1(2) I < 0,0001wr, r = ^x 2 + y 2 (11) 

and 

|Vi/i (2) | < 0,0001u (12) 
with small v (to be specified later) and these estimates are valid in the area of 
interest 

N = {y > ^} n {y < e 2 } n {x > a} 
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where 

v <C ei <C e 2 

The functions vi(2) are infinitely smooth in all variables in H but we have no control 
over higher derivatives. We also assume that the flow given by (fTOj) is area preserv- 
ing. Our goal is to study the behavior of trajectories within the time interval [0, T]. 
In this section, the parameters will eventually be chosen in the following order 

T — > e 2 — > ei — > v 

Here are some obvious observations: 
1. If £i(2) are small and 



a > v 



W (13) 
then x(t) increases and y(t) decreases. This monotonicity persists as long as the 
trajectory stays within K. Assuming that eu2) are fixed, (fl"3|) can always be satisfied 
by taking v small enough, i.e., 



v < (14) 

£2 



2. We have estimates 



x(-\ogy+C)+vy > x > x(-logy-C)-vy, -y(\\ogy\+C) < y < -y(\ logy|-C) 

(15) 

The second estimate yields 

e e\\o & P + C) >y ( t ) >e e t Qo e0 -C) (lg) 

Let us introduce 



k(T,0) 



e 



e 1 (log/3-C) 



log y(s)ds I dr 



For x(t), we have 

x(t) < a exp ^Ct — J log y (r)c?r^ + v J y(r) exp K7(t — r) — J 

x(T) < (a + vf3T) exp (t(C + | log K (T,f3) |)) 
Thus, the trajectory will stay inside H for any t £ [0, T] as long as 

a < k 3+t - ve 2 T 

and if we have 

v< K i+T (T,f3) (17) 

then the condition 

a < p Se2T (18) 
is sufficient for the trajectory to stay inside H for t € [0, T]. Thus, we are taking 



2T 



ei < e 2 

and focus on the nonempty domain 

n = {(a,P) :ei < a < f3 8e2T , f3 < e 2 } 

The condition on v is ([17| . so taking the smallest possible k(T,/3) within SIq we 
get, e.g., 

v < ef (19) 
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(20) 



Then, any point from £Iq stays inside H over [0, T], x(t) grows monotonically and 
y(t) monotonically decays with the double-exponential rate given in (|f 6p . 

Now, we will prove that the derivative in a of x(a, ft, t) grows with the double- 
exponential rate and this will be the key calculation. For any t £ [0, T], (|8]) yields 

x a = — 0.5XQ, log(x 2 + y 2 ) + x a ri + 

+xx a ri x + xy a r ly + v lx x a + v ly y a - y a arctan^y -1 ) 
y a = 0.5y a log(x 2 + y 2 ) + y a r 2 + yx a r 2x + 

+yy a r2y + v 2x x a + v 2y y a + x a arctan(yx~ 1 ) 

and x a {a, /?, 0) = f , y a {a, f3, 0) = 0. Let 

fu(t) = v Xx - 0.51og(x 2 +y 2 ) +r a +xr lx 
/i 2 (i) = xriy + v ly - arctan(xy _1 ) 
/2i(*) = yr%x + »2x + arctan(yx _1 ) 
f 22 {t) = 0.51og(x 2 + y 2 ) + r 2 + yr 2y + v 2y 



exp i^j /n(r)drj x, y a = exp 

9 = hi - f22 



t 

/ 22 (T)dr I y 



(i 



If 

then ^ ^ 

x(t) = l + J x(s)f2i{s)J fi2(r)exp(-^ g(0dA drds 

Since the trajectory is inside H, we have y > \fx and so 

|/i 2 |<W + «, l/ 2 i|<l, fn >e*(- log/3 + C), s(t)>l 
From (ITCjl . we get 

\x{ t )~ 1 \^' v j Q \x(r)\dT + J^ ^( s )l(^ e eT(log/3+c) e- (T - s) dr^ ds 
The following estimate is obvious 

" t ee ^(log^+C) e -(r- S ) dr<<e -s 



as j3 is small. Assuming that 

u<(T + l)- 1 (21) 

and e 2 is small, we have 

x(t) ~ 1 

and 

lX (e y -l)/2 



x Q (a,/3,T)> ^-j (22) 

The estimate (|22|) is the key estimate that will guarantee the necessary growth. 

Now, let us place a circle S 7 (x, y) with radius 7 and center at (x, y) into the zone 
f2o- Consider also the line segment I = [A\, A 2 ], A\ = (x— 7/2, y), A 2 = (1+7/2, y) 
in the center, parallel to OX. We will track the evolution of this disc and this line 
segment under the flow. We have by the Lagrange formula 

x(A 2 ,T) - x{A 1: T) >(3- {eT ' 1)/2 \A 2 -A 1 \ 
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From the positivity of x a (a, (3, T) it follows that the image of I under the flow is a 
curve given by the graph of a smooth function T(x). Thus, the image of I (call it 
I') has length at least (3~( e ~ 1 ^ 2 \A 2 — A\\. Denote the distance from V to S'(x, y), 
the image of the circle, by d. Then, the domain {T(x) — d < y(x) < T(x) + d, x € 
(x(A 1 ,T),x(A2,T))} is inside B' (x, y). The area of this domain is at least 

d .f}-^ T -W\A 2 -A x \ 
Thus, assuming that the flow preserves the area, we have 

Consequently, if we place a bump in D,q such that the I and Sy(x, y) correspond to 
level sets, say, h 2 and h\ (and, what is crucial, /ii(2) are essentially arbitrary 

< hi < hi < 0,0001), then the original slope of at least ~ \h% — ^i|/7 w iU become 
not less than 

r (eT - 1)/2 - (\h 2 -h x \h) (23) 
thus leading to double-exponential growth of arbitrarily large gradients. 

Remark 1. If p is a fixed small number, we have growth in T. If T is any 
positive fixed moment of time, we have the growth if f3 — > 0. 

Remark 2. Let us reiterate the order in which the parameters are chosen: we 

o 2T 

first fix any T, then small e 2l then e\ < e 2 ■ How small e 2 must be taken will be 
determined by how large the parameter A is chosen in the theorem 17.11 below. This 
defines the set Slo- For the whole argument to work we need to collect all conditions 
on v. (dl||, (HU, fl2U which leads to 

v < e{° (24) 

5. Small perturbations of a singular cross can also generate double 

exponential contraction in h 

Assume that the function 6\ at any given time t € [0, T] is such that 

0i(x,y,t) = 9 s (x,y) 

outside the "cross" -domain A = {\x — irk\ < t} U {\y — nl\ < r} where r is small 
and k,l G Z. Inside the domain A we only assume that 0\ is bounded by one 
in absolute value, is even, and has zero average. Notice here that the Euler flow 
preserves property of the function to be even. Given fixed £1(2) & n d the domain H 
defined by these constants, we are going to show that the flow generated by B\ can 
be represented in H in the form (TlOl) with v(t) — > as r — > 0. We assume of course 
that t £1 . 

For that, we only need to study 

Fi = VA- x p, p = 9 1 -6 s 

Here are some obvious properties of F\ 

1. i*i(0) = as Q\ and 9q are both even. 

2. We have 

Using the formula 

\x - y\ 



x y 

WW 
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we get 

,r|logr| 



\Fi(z)\ < \z\ 



Thus, by taking r small, we can satisfy (fTT|) . How about (JT2j ? For the Hessian, we 
have 

|JTA -1 p| < er 2 r 
and after combining we must have 

e^log^ef (25) 

by ([2"1|) . Thus, this condition on the size of the cross guarantees that the arguments 
in the previous section work. 



6. The flow generated by a small steep bump in N 

In this section, we assume that at a given moment t € [0,T], we have a smooth 
even function b(x,y,t) with support in H U — H, with zero average, and 

||&||a<w, ||V6|U<M 

(here one should think about small u> and large M). We will study the flow gener- 
ated by this function. Let 

F 2 = VA _1 6 

Here are some properties of F% 

1. F 2 (0) = 0. 

2. To estimate the Hessian of A _1 6, consider the second order derivatives. For 
example, 

(Q-e) 2 +(/3-<?) 2 <i 



l>(a-0 2 + (/3-r,) 2 >p 2 ((" - 2 V) 2 ) 



/ 



^ /Im, ft *) + V&(£', ?/, f) • (f - a, „ - /?)] 



((a-0 2 + (/?-^) 2 ) 

(q-5) 2 + (^-»;) 2 <P 2 

The first term is controlled by By our assumption, the second term is dom- 

inated by M p. Optimizing in p we have 

To guarantee the conditions that lead to double exponential growth with arbitrary 
a priory given M, we want to make u> so small that conditions (fTTj) and (fT2|) are 
satisfied with v as small as we need (i.e., (f24| ). The condition (fl2|) is immediate 
and (fTTj) follows from i<2(0) = 0, Lagrange formula and the estimate on the Hessian. 
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7. One stability result and the proof of the main theorem 

It is well known that given 0o G L°°(T 2 ), the weak solution exists and the flow 
can be defined by the homeomorphic maps ^Sg (x,y,t) for all t so that 9(x,y,t) — 
9q(^J {x, y, t)) where ^>g itself depends on 6>o- The continuity of this map though 
is rather poor ([3], theorem 2.3, p. 99). In this section, we will need to take smooth 
9q such that 

To this end, we will consider #0 = ^0 outside the domain T> of small area. Inside 
this domain we assume 8q to be bounded by some universal constant. The proof 
of Yudovich theorem (see, e.g., the argument on pp. 313-318, proof of Proposition 
8.2, 0) implies 

max maxl^OM) ~ ^e (z,t)\ -> (26) 

as \T>\ -> 0. 

This is the only stability result with respect to initial data that we are going to 
need in the argument below. 

Theorem 7.1. For any large A and any T > 0, we can find smooth initial data 9o 
so that || 0o || 00 < 2 and 

max ||V0(-,£)IU > A eT - 1 ||V0 o ||oo 
te[o,T] 

Proof. Fix any T > and find £1(2)- For larger A, we have to take smaller 62 (see 
remark 1 in the fourth section). Identify the domain J7o and place a bump (call 
it b(z)) in fio U — f2o so that the resulting function is even. Make sure that this 
bump has zero average, height ft 2 and diameter of support h\ so that the gradient 
initially is of the size ~ hijh\. Here h\ <C ft 2 <C 1 will be adjusted later. 
Take a smooth even function u)(x, y) supported on -Bi(O) such that 

u>(x, y)dxdy = 1 

T 2 

For positive small a, consider 

9 a (x, y) — 9q * uj a E C°°, Lu a = <j- 2 uj(x/<t, y/a) 

We take a <C £i so 0<t(x, y) and 0q(x, y) coincide in H. 
As the initial data for Euler dynamics we take a sum 

9„(z)+b(z) 

Then, since 0q is stationary under the flow, the stability result (|26p guarantees that 
given any r and keeping the same value of hijh\, we can find a and h\ so small 
that over the time interval [0, T] we satisfy 



C 



1. The "evolved bump" b(z, t) stays in the domain H (e.g., (t) (supp b(z) 
»)■ 

2. Outside the cross of size r (the one considered in section five) and the support 
of the evolved bump &, the solution is identical to 0g. 

Fix a and h[ so small that for any hi < h[ we have the size of A being small, i.e. 
r is small as we wish. The value of r must be small enough to ensure the double 
exponential scenario, the conditions (fTT|) and (fT2j) . For that, we need (|25]l . 
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Next, we proceed by contradiction. Assume that for all t € [0,T] we have 
||V#(2, i)||oo < M = (hi/hi)X e _1 . Then, because ||&(2, t)||2 is constant in time as 
the flow is area-preserving and ||6(z,i)||2 < hih 2 , we only need to take h 2 so small 
that \JMh\hi is small enough to guarantee the double exponential scenario and 
the estimate (|23p . This gives us a contradiction as the double exponential scenario 
makes the gradient's norm more than M (provided that e 2 <C A~ 2 ). For the initial 
value, 

II Wolloo ~ cr- 1 + hi/ hi ~ h 2 /hi 
by arranging hu2) (and keeping hi < h'j). 

Here is an order in which parameters are chosen in this construction: 

{T, A} — > e 2 — > ei — > {a, h 1(2) } 

□ 

8. The operator £ f is linearly unbounded. 

The theorem 17. f l is equivalent to the following 

Proposition 1. The operator £ t is linearly unbounded for any t > 0, i.e. 

II V£ t ^o|| 00 
sup = +00 

eoec°°(T2),o<||0o||«><iAi-U II v ^o||oo 
Proof. The proof is immediate. Indeed, given any fixed t, we have 

( ||V£ r o ||oo\ 
sup sup — — = +00 

r6[o,t] ye ec?'»(T 2 ),||9olU=i,eo-Li II V» ||oo J 

by taking A — > 00 in the theorem 17.11 Then, to have the statement at time t, we 
only need to multiply #0 by a suitable number and use remark 2 from the first 
section. 

Vice versa, in the theorem 17.11 the combination A e _1 can be replaced by arbi- 
trary large number. In this formulation, the statement follows from the proposition. 

□ 

As the statement of the theorem 17.11 holds with any A, the double exponential 
function is not relevant at all in the formulation itself. However, it is very special 
hyperbolic scenario with double exponential rate of contraction that ultimately 
provided the superlinear dependence on the initial data. 

The interesting and important question is whether the vorticity gradient can 
grow in the same double exponential rate starting with initial value ~ 1? We do 
not know the answer to this question yet and the best known bound is (see, e.g., 

0) 

max ||V0(-,t)ll°o > e aoolT 
te[o,T] 

for arbitrary T and for T-dependent 8 with ||#o||oo ~ ||V#o||oo ~ 1- 
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